Non-abelian superconducting pumps 
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Cooper pair pumping is a coherent process. We derive a general expression for the adiabatic 
pumped charge in superconducting nanocircuits in the presence of level degeneracy and relate it 
to non-Abelian holonomies of Wilczek and Zee. We discuss an experimental system where the 
non-Abelian structure of the adiabatic evolution manifests in the pumped charge. 
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If the Hamiltonian of a quantum system depends adi- 
abatically and cyclically on time via some external pa- 
rameters the wavefunction, for the system initially pre- 
pared in an energy eigenstate, after a cycle returns back 
to its initial state up to a phase, that in addition to the 
dynamical contribution, has a component of geometric 
nature [1] . It depends only on the shape of the path cov- 
ered in the parameters space and on the structure of the 
Hilbert space of the quantum system. If the adiabatic 
evolution takes place in a degenerate energy eigenspace 
the cyclic evolution leads to a superposition of the degen- 
erate eigenstates and the geometric transformation ac- 
quires a non-Abelian structure [2] . Simon [3] and Wilczek 
and Zee [2] showed that this mapping plays the role of the 
holonomy of the gauge theories. Holonomies naturally 
emerge in the description of the dynamics of simple quan- 
tum systems due to the parallel transport dictated by the 
Schroedinger equation [4, 5] which constraints the over- 
lap between the wavefunctions of the system at successive 
times to be real and positive. The holonomy group is the 
group of the transformations generated by the parallel 
transport along closed paths on the parameter space. 

Geometric effects appear naturally in adiabatic quan- 
tum pumping: in a mesoscopic conductor a dc charge 
current can be obtained, in the absence of applied volt- 
ages, by cycling in time two parameters which character- 
ize the system [6]. In the scattering approach to trans- 
port the pumped charge per cycle can be expressed in 
terms of derivatives of the scattering amplitudes with 
respect to the pumping parameters [7], the Brouwer for- 
mula [8]. Its relation to geometric phases has been elu- 
cidated in Refs.[9, 10]. If only superconducting leads 
are present pumping is due to the adiabatic transport 
of Cooper pairs. Besides the dependence of the pumped 
charge on the cycle, in this case there is a dependence on 
the superconducting phase difference(s) (the overall pro- 
cess is coherent). Cooper pair pumping has been studied 
both in the limit of transparent interfaces [11, 12] and 
in the Coulomb Blockade regime [13-15]. A connection 



between Berry phase and pumped charge has been es- 
tablished also in this case [11, 16, 17] thus opening the 
possibility to detect geometric phases in superconduct- 
ing circuits [17, 18]. An experiment of this kind has been 
successfully performed recently in Ref. [19] thus paving 
the way to holonomic quantum computation [20, 21] with 
superconducting nanodevices. 

In this Letter we study Cooper pair pumping in super- 
conducting circuits in the regime of Coulomb blockade. 
The new feature we consider here is the possibility to 
pump in a degenerate subspace. We derive an expression 
for the pumped charge in the presence of a degenerate 
spectrum and relate it to the non-Abelian connection of 
Wilczek and Zee. We propose a superconducting net- 
work where this relation can be tested and discuss two 
clear signatures of non-abelian holonomies. First, under 
appropriate conditions, the pumped charge per cycle is 
quantized. Second the pumped charge depends both on 
the cycle and on the point where the cycle starts. If 
tested experimentally this would be a clear proof of the 
non-Abelian nature of pumping. 

The possibility to generate non-abelian holonomies in su- 
perconducting circuits has been studied previously [22]. 
Ref. [22] dealt with the problem of implementing a holo- 
nomic quantum computer with Josephson circuits. In 
that work the authors show that a cyclic adiabatic change 
of parameters in a closed system may lead to a final state 
having a different charge distribution than the initial one. 
In that context the word "charge pumping" was used to 
describe such process even if the charge cannot be re- 
ally pumped in or out the system. In the present work 
we study charge transport between two external reser- 
voirs connected to the circuit via Josephson junctions, 
the word "pumping" has thus a different meaning. In 
order to keep the presentation transparent we discuss 
an idealized situation, in the conclusions we discuss the 
various problems that may occur in experiments. The 
Cooper pair pump consists of a Josephson network con- 
nected through Josephson junctions to two superconduct- 
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ing leads. An example (at this stage the discussion is 
general) is presented in Fig. 1. The two superconduct- 
ing electrodes are kept at a finite phase difference = 
Vl ^ where '^rt,/L is the phase of the superconduct- 
ing order parameter of the right/left lead. The Cooper 
pair pump is operated by changing adiabatically in time 
some external parameters such as gate voltages or mag- 
netic fluxes. We will label this set of external parameters 
by the vector \{t) = {Vgi{t),^i{t)}. The Hamiltonian 
of the pump depends on the superconducting phase of 
each island of the network (pj (z = 1, . . . A^), on its conju- 
gate, the charge on each island n^, on the phase difference 
across the pump, and on all the external parameters: 
H{i) = H\^pl, ...,(fN',ni,...,nN] A(t), ip] . The state of the 
system is denoted by |\I/(t)) and |\E'(t)) = A(i), tp] ). 
By changing the control parameters in time, a charge 
qC*'') will be transferred in the circuit. The transferred 
charge after a period T can be obtained by integrating 
the charge imbalance between the outer capacitors that 
connect the network to the leads 

Q(tr) ^ _2ig r Q^, [^^(^/-,| |^(^/)^] _ (1) 

Jo 

This definition of the transferred charge may be derived 
from the time integral of the current operator, Q^*''^ = 
{^{t')\{d^H{t))\^{t'))dt', and the Schrodinger 
equation. In the definition of the transferred charge 
is not a quantum operator but a simple derivative re- 
spect to a classical parameter (the phase difference be- 
tween the two electrodes). In this Letter we generalize 
the results obtained so far relating Cooper pair pump- 
ing to geometric phases allowing the spectrum of II{t) to 



be degenerate. We assume that for all A in the parame- 
ter space a degenerate energy eigenspace Tin (A) exists of 
constant dimension, corresponding to the eigenvalue 
En (a) . The control parameters are varied in time adia- 
batically, i.e. A = X{t/T). It is convenient to introduce 
for all t e [0,T] a basis, {ipna{t)} a G [l..D„], of the de- 
generate subspace, 7i„(A), formed by the instantaneous 
eigenstates of the hamiltonian, H{t): H(t) |V'fea(^)) = 
Ek{t) li'kait))- As discussed by Wilczek and Zee [2], if 
initially the state of the system is in one of the degen- 
erate eigenstates = 0)) = |'0„^[A(O)]> e W„(A(0)), 
then after a cyclic evolution 

|*(T)) = [t/„(T)]„, |Vna(A(0))) + 0{l/T) . (2) 

The Dn X Dn operator Un{t) can be written as 

Un{T) = e-*r^"WTe[-r'^"M . (3) 

Here, T denotes the time-ordering and the connection 
r„(i) is given by [r„(f)]^^ = {tpna{t)\ ipnoit)^- The rela- 
tion between the transferred charge and the non-Abelian 
holonomy can be obtained by substituting (2) in (1) 
[23]. We assume that the system is prepared at f = 
in a linear superposition of the degenerate eigenstates, 
I ^ira) = C7I V'n7[A(0)]); duriug a cycle of duration T 
the total transferred charge then reads 

Q<*'^^=E<^7'Q77" (4) 
77' 

where the charge matrix / is given by 
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The first term in the r.h.s. is the supercurrent contri- 
bution to the transferred charge. The second and third 
terms are of geometrical nature and describe pumping. 
If the pumping occurs through a non-degenerate level, 
Dn = 1, we recover the abelian result [16, 17, 24]. The 
measurement of non abelian holonomies is a non-trivial 
task, in particular the non commutativity of the theory 
has revealed not easy to detect [5] . Non-abelian contribu- 
tions to pumping introduce qualitatively new effects that 
can be verified experimentally. We will highlight these 
aspects by analyzing a specific case that can be realized 
experimentally. 

A possible experimental realization of a non-abelian 
pump is shown in Fig.l. It is a three island pump with 



four symmetric Josephson SQUID loops: two inner loops 
with capacitances Ci and C2 and two outer loops with ca- 
pacitances Cl and Cr. The outer loops connect the net- 
work to the superconducting electrodes which are kept at 
a constant phase difference <p (from now on for simplicity 
we fix (^i = and <^_r = 0). The charging configuration 
of the system can be controlled externally modulating 
three gate voltages, Vu, Vgi and Vg2 connected to the is- 
lands via the respective gate capacitances, C„, Cgi and 
Cg2- The effective Josephson couplings, J^, Jr, Ji, J2, 
can be tuned independently varying the magnetic fluxes 
through each loop. All the Josephson coupling energies 
are much smaller than the charging energy of the sys- 
tem, Ec- The charge states are indicated as |n„,ni,n2). 
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rii being the excess charge on the i-th island in units 
of 2e. The reahzation of the degenerate subspace re- 
quires additional constraints on the gate voltages. We 
take the gate capacitances to be small compared to the 
Josephson capacitances, i.e. Cu ~ Cgi <^ Ci ^ Cj with 
i — 1,2 and I ~ L,R, and we assume that the gate charge 
rigu = CuVgu/2e is in the range i < ng„ < | + Cl+Cr 
while the two gate charge Ugi = CgiVgi/2e satisfy the 

condition Ugi = i ^1 + ^(2ng„ — 1)^ with i = 1, 2 and 

Ct = Ci + C2 + Cl + Cji- Under these conditions 
only four charge states are relevant to the dynamics at 
low temperatures: |0,0,0}, |0, 1, 0) , |0, 0, 1) , 1 1,0,0) (ah 
other charge states are at a much higher energy Ec)- 
Moreover the charge states |0, 0, 0} , |0, 1, 0} , |0, 0, 1) are 
degenerate while the charge state |1, 0, 0), corresponding 
to the configuration in which there is one excess charge 
on the island U, has higher energy. In this restricted 
Hilbert space the effective hamiltonian of the pump can 
be written as 

H = S„|l,0,0)(l,0,0| + [Jeff(</p)|l,0,0)(0,0,0| 
+ Ji|l,0,0}(0,l,0| + J2|l,0,0}(0,0,l|+h.c.] (6) 

where Jctfif) = {Jl^-^"^ + Jr)- We set to zero the 
electrostatic energy of the degenerate charge configu- 
rations and we denoted with Eu the charging energy 
of the state |1,0,0), E^ = ^''(g^+pffi^^ ■ Given the 
capacitances, our choice of the gate voltages guarantees 
that the charging hamiltonian is symmetric under 
the simultaneous exchange of the three charge states 
and of the three couplings Jeff, Ji and J2. This fact 
leads to a two dimensional degenerate subspace for 
any value of the couplings. Hamiltonian (6) has been 
discussed previously in the context of adiabatic passage 
techniques in a quantum optics [25] and supercon- 
ducting nanocircuits [27] and for holonomic quantum 
computation as demonstrated in [26]. The crucial point 
here is that Eq. (6) arises from a Josephson network 
in the presence of superconducting electrodes. This 
is why we are able to relate pumping to holonomies. 
The hamiltonian has three distinct eigenvalues: £^0 — 
1/2 [E^ - {El + 4( J2 + Jl + J|)) 1/2] ^E^ = Q and E^ = 
1/2 [Eu + [El + 4(J2 + j2 + j2))i/2] . In the previous 

1 /2 

definitions we set Jq — {Jl + Jr + "^JlJr cos ip) . The 
eigenvalue Ei remains zero and it is doubly degenerate 
for any value of the Josephson couplings not all zero. An 
orthonormal basis for the two-dimensional degenerate 
subspace corresponding to the eigenvalue Ei is given 
by : IV'u) = iVii ( J2 [0, 1, 0) - Ji |0, 0, 1}) and IV'12) = 
^^12 [{J! + Ji) |0, 0, 0) - Jeff ( Jl |0, 1, 0} + J2 |0, 0, 1))] 
where Nn and N12 are normalization factors. 

In the following we will show that: (i) The pumped 
charge may be quantized; (ii) By composing two cycles in 
the parameter space, the pumped charge depends on the 
order in which the two cycles are followed, (iii) Consider- 
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Figure 1: Non-abelian superconducting pump. The gate volt- 
ages are kept fixed during the cycle, their values are chosen 
in order to have a doubly degenerate spectrum. The only 
pumping parameters are the fiuxes through the four loops. 



ing the pumped charge as a function of time, Q = Q{t), 
the period of Q{t) may be a multiple of the period of the 
parameters cycle. The second and the third point are due 
to the non-abelian character of the adiabatic evolution. 

(i) We first discuss the quantization of the pumped 
charge. To this end the cycle can be divided in three 
steps. First, we prepare the system in the charge state 
0,0,0) and we set all the couplings to zero except for 
J2 which is kept constant and positive during the whole 
cycle. For these initial values of the couplings the 
two degenerate eigenstates are also charge eigenstates: 
\^ii{t - 0)) = |0, 1, 0), |V^i2(t = 0)) = |0, 0, 0); so the sys- 
tem is initially in the eigenstate l'0i2). Second, we per- 
form a 7r/2 rotation in the degenerate subspace by manip- 
ulating adiabatically and cyclically Jl and Ji and keep- 
ing Jfl to zero. In this phase one charge enters the pump 
from the left reservoir and the state vector of the system 
undergoes the transformation: [0,0,0) — > e"^ [0,1,0). 
Third, we let the charge out of the circuit with another 
7r/2 rotation in the degenerate subspace. During this 
third phase we manipulate Jr and Ji and we keep Jl 
to zero. The state of the system undergoes the transfor- 
mation e''^ [0, 1,0) e'^f [0,0,0). How is it possible to 
realize the 7r/2 rotations? One can show that by manip- 
ulating adiabatically Ji,J2 and Jl or Ji,J2 and Jfi and 




Figure 2: By varying the Josephson couplings through the 
cycles (^L (left) and then (right) one Cooper pair is trans- 
ferred through the pump. 
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1,(0=0 



by: Q*'- 



(3^ 



Ji(t)=J 

jR(t)=0 



Figure 3: Non commuting cycles in the parameters space. 
Along ^1, J2 is set to zero while along (^2, both and Jr are 
fixed. The circles indicate the initial points of each cycle. 



keeping respectively Jr or J^, zero the holonomy reads 



r/(C) 



^VL/R sin 



.g-'Vi/Rsinf]^ cosf2^ 



(7) 



where Sl,^ = /^(V'li |V'i2)- The required rotations are re- 
ahzed whenever fi^ = 7r/2, for example by means of the 
cycles Cl and Ct?, shown in Fig. 2. By substituting the 
evolution U , the connection F and the matrix elements 
of {i}na\d,p\ipni3) for the cycles of Fig. 2 in Eq. (5) we find 
Q*^*''-' = — 2e. For the present cycle, the charge pumped 
starting from the eigenstate IV'ii) has opposite sign to the 
charge pumped starting from the eigenstate |V'i2)- (ii) 
The non-abclian nature of the evolution has measurable 
effects on the pumped charge. We consider the cycles de- 
picted in Fig. 3 and determine the pumped charge in the 
two cases in which the pumping cycle is cither performed 
by first going over ^1 and then ^2 or in the reverse order. 
In the abelian case considered in Refs. [14, 16, 17] the two 
situations arc equivalent. In our case the order matters, 

i.e. there are examples where Q^*'^^ 7^ ^^III-l ■ ^ specific 
example is provided by the pumping cycle obtained by 
performing the two loops presented in Fig. 2. Assum- 
ing that the initial state is the state jO, 1, 0), the pumped 
charge in the two cases is Q^q^q^ = § and Q^^^^^ = e. The 
pumped charge differs in the two cases. 

(iii) Another manifestation of the non- Abelian nature 
of pumping is that after a cyclic evolution of the external 
parameters the state docs not necessarily go to the ini- 
tial state (sec Eq. (2)). One can therefore design paths 
in parameter space such that after N cycles the system 
returns, up to a phase, to its initial state. In this situ- 
ation the pumped charge per cycle will not be constant 
in each cycle but it will have a period which is NT, T 
being the period of the elementary cycle. In fact, as- 
sume that a certain cycle is performed consecutively 
N times starting with the system in the state \^in) de- 
fined before. At the beginning of the N-th cycle the state 



will be (C^^)^ ^ l^'iri) then, as one can easily show us- 
ing Eq. (4), the total transferred charge will be given 



Nth cycle 



= E 



77' 



JV-l 



Eventually, when ^ is proportional to the identity 
operator, the transferred charge will be periodic with pe- 
riod NT. 

The situation we discussed so far is ideal. Several im- 
portant issues have to be considered in a realistic situa- 
tion. First of all, it would be desirable to pump through 
the ground state; here the degenerate subspace is the 
first excited level, which is sensitive to relaxation effects. 
We do not think this is a serious problem, though: we 
showed how the state can be accessed and moreover the 
real bottleneck is the decoherence time and not the re- 
laxation time, which typically is much longer. Indeed, 
the effect of decoherence on coherent pumping appears 
to be a more fundamental issue since the presence of an 
external bath may, in addition to relaxation, lift part of 
the degeneracy which is crucial for non- Abelian pump- 
ing. In addition, the degeneracy may be lifted because of 
the unavoidable static imperfections in the network. We 
do not expect degeneracy lifting to prevent the observa- 
tion of non- Abelian effects on pumping, it just imposes 
a constraint on the duration of the cycle: T should be 
shorter than min{?i/AiJ, Ti/r^} where /S.E is the small 
splitting arising from the non perfect degeneracy of the 
levels involved and the decoherence time. 
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